Spin Josephson vortices in two tunnel coupled spinor Bose gases 
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We study topological excitations in spin-1 Bose-Einstein condensates trapped in an elongated 
double- well optical potential. This system hosts a new topological defect, the spin Josephson vortex 
(SJV), which forms due to the competition between the inter- well atomic tunneling and short-range 
ferromagnetic two-body interaction. We identify the spin structure and formation dynamics of the 
SJV and determine the phase diagram of the system. By exploiting the intrinsic stability of the 
SJV, we propose a dynamical method to create SJVs under realistic experimental conditions. 
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Ultracold spinor atomic gases that exhibit both super- 
fluidity and magnetic order display an abundance of rich 
static and dynamical properties. This has attracted con- 
siderable theoretical and experimental study, with par- 
ticular focus on the topological excitations of trapped 
spinor gases [THS]. Topological phases of single trapped 
spinor gases, such as spin vortices [6H9], knots [10] and 
skyrmions [TTHT3] . depend critically on the mean- field 
order-parameter manifold. A remarkable feature of these 
dynamical excitations is that their size is typically larger 
than the underlying spin healing length (SHL) [TJ[2]. 

When spinor atoms are confined in optical lattice po- 
tentials, atomic tunneling between adjacent lattice sites 
competes with the spin dependent inter-atomic interac- 
tion. This competition provides a mechanism for the 
emergence of topological phases, which have been iden- 
tified and investigated in several studies [T4H2T] . It can 
also strongly influence the behavior of a simpler system 
comprising atoms confined in double-well (DW) poten- 
tials [22j|23], which are analogous to Josephson junctions 
in solid state devices. Analysis of such systems often uses 
the lowest energy mode approximation [24] , which allows 
the intra-well spatial motion to be mapped as a function 
of time. Even in this limit, spin-dependent population 
oscillations between the two potential wells have been 
identified [25H29] . But beyond this limit, it remains un- 
clear whether quantum fluctuations can trigger the for- 
mation of extended topological excitations when the size 
of the individual spinor gases in each well exceeds the 
spin healing length (SHL). 

In this work, we show that a dynamically stable topo- 
logical excitation, the so-called spin Josephson vortex 
(SJV), forms in two weakly coupled spin-1 ferromagnetic 
Bose-Einstein condensates (BECs) trapped in an elon- 
gated DW potential [Fig.JTJa)]. As depicted in Fig. [TJb), 
a key feature of an SJV is its fixed spin current, facilitated 
by the inter- well atomic tunneling, which circulates about 
a point mid-way between the two wells. Due to its large 
size, on the order of several SHLs, the SJV is a macro- 
scopic topological object. We determine analytically the 
parameter space required for SJVs to form in a uniform 
system, where they are the only stable topological exci- 
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FIG. 1: (Color online) (a) Schematic diagram of the system. 
The weakly coupled spin-1 BECs are trapped in a double- 
well potential, which gives strong (weak) confinement in the 
y — z plane (x-direction). A uniform magnetic field Bz is 
applied along the z-direction. (b,c) show spin vector patterns 
corresponding to (b) a spin Josephson vortex (SJV) centered 
on the blue cross, (c) a ferromagnetic domain wall (FDW). 
The parameters are k, = k, c /2 for the SJV and k, = 2k c for the 
FDW. Other parameters are q = and a = tt/2. See text for 
more details of the parameters and spin patterns. 



tation. We show that, as a consequence of this stability, 
the SJV can be created dynamically through the decay 
of a ferromagnetic domain wall (FDW) [6] [30j [31] . We 
demonstrate that the SJV can be realized by implement- 
ing this dynamical scheme under conditions that can be 
fully attained with current experimental techniques. 

Our system comprises two one-dimensional (ID) spin- 
1 BECs trapped in a symmetric optical DW potential 
[Fig. [TJa)]. Both atom clouds are strongly confined in 
the transverse (y, z) directions. The dynamics of weakly 
coupled spinor BECs may be described by the spin- 
1 field operator 32]. 4f(x) = 4? l (x) + 4? r (x), where 
4?i(x) = [V^i(x), ^o( x )> ^-i(x)] T with j = I (r) indicat- 
ing the left (right) well and m = {1,0, —1} denoting the 
three Zeeman levels. A uniform magnetic field Bz is ap- 
plied along the z-axis. The many-body Hamiltonian is 



H = H t + H l +H r , 



(1) 



where % t = — « / dx\^f l ^ r + H.c] denotes the inter- well 
tunneling with state-independent tunneling strength n 



2 



[24] . The Hamiltonian of atoms in the j th well is W = 
/da;[*^W'^' + f (n^) 2 + fF^'-F^, where h^ m , = 
-£ m , m '(<9 2 /2 - fi j - prn + gra 2 ). Here, fi j , p = -/jLbB/2 
and g = (ii&B) 2 / AE^ denote the chemical potential, lin- 
ear and quadratic Zeeman energies respectively, where 
/iB is the Bohr magneton and E^f is the hyperfine en- 
ergy splitting [31 j. The two-body collisional interactions 
in the j th spinor BEC enter the expression for H? via 
the scalar density n J = ^t^rJ anc [ the spin-dependent 
vector density, F J = ^'tf ^ro where f is the Cartesian 
vector of the spin-1 matrices (f x ,fy,fz)- The effective 
ID interaction strengths are cq = 16ft 2 (ao + 2a2)/9Mr 2 L 
and c\ — — 16ft 2 (ao — a2)/9Mr 2 L , where as is the 3D s- 
wave scattering length for collisions with total angular 
momentum S = 0,2 [TJ [2], and r± is the width of the 
BEC in the transverse directions. 

We study the system using mean field theory. Let us 
first investigate the stationary state of the system. We 
use a simple ansatz to describe the order parameter 
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where n J = ^t^ri i s the total density of the j 
BEC [31 j. For convenience, we scale length, time and 
energy by the variables £o = ^/(Mcon^) 1 / 2 , to = h/c^riR 
and eo = c^ur, where the reference density, is 
chosen to ensure correct chemical potentials in the two 
atom clouds. Neglecting the spatial dependence of n J , 
if){ satisfies the coupled non-linear differential equations, 
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0, where 7 = ci/cq 

and j = I (r) when j' = r (I). Consequently, within 
this approximation, the spinor BECs are described by 
two coupled scalar equations. Each scalar equation is 
characterized by an effective chemical potential /x; 
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(1 + 27) n 3 — q and an interaction strength equal 
to —47 > 0. The corresponding stationary solution is 
readily obtained [33 



V>i = [Ctanh(ra) ± iAsech(vx)} e % 



(3) 



where C = \J (2yn) + q)/8 / y, v and A are constants and 
the +(— ) sign corresponds to j = / (r). This analytical 
solution allows us to calculate many properties of the 
system. For example, one can directly find the density 
n J = (1 + ac + 7)/(1+7) in units of ur and chemical 
potential \J = q/2 + (7 + 1) in units of eo- 

We now discuss the topological excitations of the sys- 
tem. Depending on the value of A, two distinct solutions 
can be obtained from Eq. When A = 0, the solution 
describes a FDW [31 , which has a characteristic spatial 
width of \ jv = V2£, s , where £ s = l/(2|7|n J — g)s 5 is the 
spin healing length. When A = yj (2^ni + q + 8k) / 87, 
we obtain a totally different topological excitation, the 
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FIG. 2: (Color online) (a) and (b): \^[{x)\ 2 for a FDW and an 
SJV respectively using the same parameters as in Fig.[l] The 
dashed (dotted) curves are calculated with (without) the con- 
stant atom density approximation and the solid curves show 
numerical solutions found by evolving the GPEs in imaginary 
time, (c) Phase diagram of the topological excitations. The 
SJV is dynamically stable in region I (yellow). The FDW is 
unstable in both regions I and II. In region III, the system ex- 
hibits a polar groundstate phase (see text). Red (blue) circle 
marks system parameters in region I (II) , which are discussed 
in the text. 



spin Josephson vortex. The size of an SJV is approx- 
imately 1/v = l/y/4k~ and, hence, controlled by the 
inter-well tunneling strength, n. To distinguish the 
two distinct topological excitations, we calculate their 
spin texture, characterized by the local spin orienta- 
tion (j)i(x) = tan -1 (FJ/FJ) and its magnitude |F J | = 
[(FJ) 2 + (i^') 2 ] 1/2 , from Eq. The spatial variation 
of the local spin vector along the x axis is shown in 
Figs.[jjb,c). In an SJV [Fig.Qb)], the spin current forms 
a vortex structure in which the local spin vector rotates 
between the two spinor BECs around a point [blue cross 
m Fig. [Tpb)] mid -way between them. By contrast, there 
is no spin current associated with the FDW. Instead, the 
spin vectors in the two atom clouds are locally aligned 
for all x and vanish at x = [Fig. [TJc)]. 

Although the analytical ansatz in Eq. Q is simple, 
it produces accurate wavefunctions when compared with 
full numerical solutions of the equations of motion, which 
we obtain by propagating the Gross-Pitaevskii equa- 
tions (GPEs) for the coupled spinor BECs in imaginary 
time [34 . Figs. [2ja,b) reveal a small deviation between 
ip[(x) curves obtained analytically (dashed curves) and 
numerically (solid curves) near the center of the SJV and 
FDW. This deviation is caused by the constant density 
assumption used in the above analytical calculation. To 
overcome this, we now allow a spatially dependent den- 



3 







U 

10 

on 
ZD 




—SO 


\tft x 
ft 1 


— 40 


50 




60 









10 




20 
^30 


tttf 


—40 
50 
60 



(6) I 

jTWT\ 


(c) ' 


(d) 1 






(e) ' 








(/) 



-2-1012 -2-1012 



FIG. 3: (Color online) (a) Schematic diagram showing phase imprinting of the coupled spinor BECs by a focused laser beam 
(red). The phase-imprinting laser is switched on at t = and only affects atoms in the region — rb < x < Tb spanned by the 
beam (upper panel). The laser beam is switched off at time r = tt/(/3Io) and coherently flips all atomic spins in the region 
— rb < x < Tb (lower panel), (b) and (c) Color maps showing how the spin vectors and F r , respectively, evolve after the phase 
imprint when k = 2k c . The orientation (<jp) and magnitude (|F J |) of the local spin vector in the x — y plane are represented, 
respectively, by the color and brightness of the images (see scale). The two black stripes (where |F J | = 0) visible for \j\t < 5 
at x/rb — ±1 show the initial formation of two FDWs. At later times the stripes vanish, indicating the decay of the FDWs. 
(d) and (e) are the same as (c) and (d), respectively, except that k = « c /2. Again, two black stripes centered at x/rb = ±1 
indicate the initial formation of two FDWs for \j\t < 5. In this case, though, the FDW evolves towards a new quasi- static spin 
texture [green and red stripes in (d) and (e)]. At \j\t = 50, this spin texture corresponds to the local spin vectors shown in 
(f). Comparing the region of (f) within the blue dashed box to Fig. 1(b), we see that an SJV with a = tt/2 (see text) has 
formed, centered at x/rb ~ —1 (blue cross). An anti-SJV with opposite spin vector rotation has formed within the red dashed 
box centered at x/rb ~ 1 (red cross). In our simulations, r x = 1250 and rb = 100 (see text). 



sity perturbation, n j (x) = n J +(5n J (x), in the ansatz. The 
resulting values of ip[(x) [dotted curves in Figs. [2|a,b)] 
agree much better with the values obtained numerically 
(solid curves). 

We are now in a position to obtain the phase diagram 
of the topological excitations from the above stationary 
solutions. Three distinct phase regions are found, which 
are summarized in Fig.[2|c). Region I (yellow) shows the 
parameter space where the SJV exists. The FDW exists 
both in region I and region II (green) and a "polar" phase 
in which all atoms are in the m = spin level occupies 
region III (blue) [TJ [35] . We emphasize that the phase 
diagram can be determined completely analytically. For 
example, the phase boundary between region I and II oc- 
curs along the red dotted line in Fig.[2|c), whose equation 
is« = K c = (2|7|- g )(l-| 7 |)/(4-6|7|). 

To investigate the dynamical stability of the topo- 
logical excitations, we use an extended Bogoliubov the- 



ory [36] in which we evolve a stationary solution, *&i(x), 
to & j (x,t) = ^i(x) + 6& j (x,t) at time t, where 
S^^x^t) = ui(x)e~ lXt — v J (x)*e* A * t is a small pertur- 
bation. Linearizing these vector equations with respect 
to u J and v J yields an eigenequation with eigenvectors 
(u z , v*, u r , v r ) T and eigenvalues A. A stable solution re- 
quires that Im(A) = [37]. Analysis of the eigenvalues 
reveals that the SJV is dynamically stable in region I. By 
contrast, the FDW is unstable in all regions of Fig. [2|c). 

Guided by this stability analysis, we now explain how 
to realize the SJV in experiment. First, a FDW is cre- 
ated in the spinor BECs by a phase-imprinting method 
of the type used previously to generate topological exci- 
tations [9] [38j [39] . Provided the system is in region I of 
Fig.^c) (k < k, c ), the unstable FDW can decay into the 
stable SJV. 

We now consider the details of the spin-dependent 
phase-imprinting process. A phase-imprinting laser beam 
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propagating along z is switched on at t = 0. We approx- 
imate the intensity profile of the beam along the x di- 
rection by a square wave of the form I(x) = Io6(x ± r^), 
where Iq is the laser intensity, is the Heavyside function 
and 2r\) determines the width of the laser beam along the 
x direction. Such a shape can be achieved, for example, 
by reflecting the laser beam from a spatial light modu- 
lator [39 . The laser light is circularly polarized (cr+) to 
induce a linear Zeeman shift through the spin-dependent 
A.C. Stark energy shift p(x) = /3I(x), where f3 is a con- 
stant [40 j. Applying the beam for a duration r = 7r/(/32o) 
coherently flips the atomic spins near the central region 
of the atom clouds where \x\ < rv Such a process is 
shown schematically in Fig. [3ja) . 

As in typical cold atom experiments, we further as- 
sume that the spinor atoms are confined along the x di- 
rection by a shallow harmonic trap. The atom density 
of the spinor BECs is then given by a Thomas-Fermi 
profile, n(x) = ur[1 — (x/r x ) 2 ], where r x is the Thomas- 
Fermi radius. Initially, the two BECs are prepared in 
the ferromagnetic groundstate, whose wavefunction is 
&> r = y/n(x)[-l/2,i/V2,l/2] when q = 0, in which 
the spin vectors point along the y direction pQ. After 
applying the phase-imprinting laser, we determine the 
dynamics by evolving the coupled GPE for the spinor 
BECs. We include dissipation in our numerical simula- 
tions following the methods in [8]. 

Let us now analyze the dynamics of the atom clouds 
after the laser illumination. We first consider the evo- 
lution of a system with parameters located in region 
II of Fig. [5|c). Specifically, we choose q/\j\ = and 
k/\j\ = 2k c /\j\ [marked by the blue circle in Fig. [2|c)], 
for which the evolution of the local spin vectors in the 
left and right atom clouds is shown in Fig. [3^b) and (c) 
respectively. For short times (|7|t < 10), the system re- 
acts to the laser by quickly forming two FDWs at the 
edges of the phase-imprinting region {x/r^ — ±1). These 
appear in Figs. [3jb,c) as dark stripes, where the mag- 
nitude of the atomic spin vectors |F J | = 0. The stripes 
separate two bright yellow regions, where the atom spins 
point along the y axis, from a bright blue region where 
the atom spins point along the —y direction. Due to their 
instability, the FDWs decay into spreading spin textures 
when \j\t > 10. This decay appears in Figs. [3^b,c) as 
muti-colored bands, which emerge from the black stripes 
and spread outwards with increasing t. 

The evolution of the spin texture differs markedly 
when the system parameters are prepared in region I 
of Fig. 2jc). To illustrate this, we choose q/\j\ = 
and ft/ 7 1 = ft c /(2|7|) [marked by the red circle in 
Fig. [2|c)], for which the evolution of the local spin vec- 
tors is shown in Figs. j3^d,e). Comparison of Figs. (3^d,e) 
with Figs. |3jb,c) shows that in both cases the system 
initially (for \j\t < 10) forms two FDWs (black stripes 
at x/rb = ±1). However, instead of decaying into a 
spreading spin texture, in Figs. [3|[d,e), the FDW evolves 



towards a different quasi-static spin pattern, which ap- 
pears as two vertical bright green [Fig.[3jd)] or bright red 
[Fig.^e)] stripes centered at x/rb = ±1. To demonstrate 
that these spin textures correspond to SJV formation, in 
Fig. [3jf) we show the associated spin vector configura- 
tions in the two spinor BECs at \j\t = 50. Comparison 
of the spin vectors around the point x/r^ = —1 (within 
the blue dashed box) with Fig. [TJb) clearly shows that an 
SJV with a = 7r/2 has formed. Rather less obviously, a 
so called anti-SJV has formed at x/r^ = 1 (within the red 
dashed box). This corresponds to a solution in Eq. ^ 
with a = — 7r/2. 

In practice, we can create the SJV using, for example, 
87 Rb BECs. If the total number of atoms is 2 x 10 6 and 
fx(f±) — 200 jam (2.4 jam), the characteristic timescale is 
|7|£o = 16 ms. The atomic tunneling strength, ft, can 
be controlled by changing the intensity and/or waist of 
the laser that creates the double- well trap [23] . All of the 
system parameters and procedures required to implement 
our proposed route to creating SJVs can be attained us- 
ing current experimental setups [35]. Consequently, we 
expect that the dynamical regime that we have identified 
will be directly accessible to experimental study. 

In conclusion, we have identified SJVs in spin-1 ferro- 
magnetic BECs trapped in an elongated DW potential. 
We have presented a detailed analysis of the stability and 
formation of the SJVs. In particular, we have shown that 
the SJV can be created from the decay of a FDW. Our 
analysis can be extended to study topological phases in 
mult i- well optical potentials and for higher atomic spins, 
which seem certain to reveal further exotic spin textures. 

This work is funded by EPSRC. WL acknowledges 
funding through an EU Marie Curie Fellowship. 
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